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ABSTRACT 

The paper describes a procedure for modelling the anisotropic elastic-plastic 
behavior of metals in plane stress state by the mechanical sub-layer model. In this 
model the stress-strain curves along the longitudinal and transverse directions 
are represented by short smooth segments which are considered as piecewise linear 
for simplicity. The model is incorporated in a finite element analysis program 
which is based on the assumed stress hybrid element and the viscoplasticity-theory. 

1 . INTRODUCTION 


For time-independent elastic-plastic behaviors a very convenient model to 
represent kinematic hardening is to use an assembledge of elastic-perfectly-plastic 
elements to represent the stress-strain relation which is approximated by a curve 
with several piecewise linear segments. This model, which has been widely used 
for numerical analysis of multiaxial elastic-plastic behavior is named mechanical 
sublayer or overlay models [refs. 3,4,5]. For more general case including three- 
dimensional solid, the method should perhaps be called mechanical sub element 
method. 

For plane stress problems the corresponding mechanical model is a laminated 
plate with layers of elastic-perfectly-plastic materials of different yield 
stresses. Differential equations for the solutions of plane stress elastic-plastic 
and isotropic stress-strain relations has been obtained for model with two layers, 
one of which is elastic and the other is elastic-perfectly-plastic [ref. 6]. The 
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equations are non-linear and for the case of uniaxial loading the resulting strain 
hardening behavior will not be a straight line. Thus, for a material with uniaxial 
stress-strain relation approximated by straight line segments, it is, strickly 
speaking, not possible to obtain exact representation by a mechanical sublayer 
described above. In Ref. 6, a relationship has been obtained between the ratio 
of the initial tangant modulus and the elastic modulus and the thickness ratio of 
the two layers. It is however, reasonable to assume that by using sufficiently 
small segments a piecewise linear model can be adopted. 

Hunsakier et al. [ref, 5] have also obtained a corresponding relationship for 
three-dimensional isotropic solids. In that case, for a material represented by one 
elastic subelement and one elastic-plastic subelement, the resulting uniaxial stress- 
strain relation will have linear strain hardening behavior. The proportion of the 
volume of the elastic-plastic element to the total volume V.j/V is expressed simply as 



where is the elastic modulus and £ 2 * the tangent modulus. In this case, when 
the uniaxial stress-strain relation is represented by linear segments a corresponding 
subelement model can be constructed exactly. 

The present paper is to extend the mechanical sublayer model to materials 
with anisotropic plastic behavior. Again the plane stress problem is considered. 
Finite element method for elastic-plastic analysis based on the viscoplasticity 
theory and the stress hybrid model is used in conjunction with the present mechanical 
subelement model. An example solution of a time-independent elastic-plastic problem 
is presented. 
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2. Mechanical Sublayer Model for Anisotropic Plasticity Problems 

Figure 1 is a plate with two layers under plane stress loading. Layer 1 is 

elastic-perfectly-plastic and is considered tranversely isotropic with yield stresses 

Y and Y respectively along the longitudinal and transverse directions. Layer 2 is 
X y 

elastic. The elastic constants E and Poisson's ratio v for both layers are identical. 
The yield condition for layer 1 is governed by the Hill's generalized yield criterion 




( 1 ) 


where a = Y under uniaxial loading along x direction. With yield stresses under 

X 

uniaxial loading along y and z direction equal to Y^, we can express the constants 
F, H etc. in terms of the yield stresses and obtain the following yield conditions 
for the plane stress problem 


f=[a^-a a + a o ^ + a o - Y = 0 

x^ y^ y^ s xy^-* x 


where 


a = (Y /Y f 
x' y^ ^ 

a = (Y /Y y 
s ' X xy' 


The flow rule is 
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Now, consider the behavior of the laminated plate under inplane loading 

conditions. From testing tension coupons cut along the x and y directions, the 

elasticmodulus for both direction is given as E-j while the initial tangent moduli 

are represented, respectively, by E and E as shown in Figure 2. Let a and o 

^2 '^2 ^ ^ 
represent the average in-plane stresses, then the stress rate is 

t, to 

where t is the total thickness and t-j and t 2 are the individual thickness. The 
thickness ratio then is. 


Here, layer 2 is elastic hence. 


(Ey 


From f = 0, we obtain. 


(2 a - a„ ) + a, (2ao„ = ) = 0 


xi y-, 


yi X., 


Thus, at initial yield when a = o =o = 0 

y y 1 y 2^ 


and from Eg. (4), 


We also have 


= 2 

“^1 


• P - o’ ^ 

e.. 


C - {o -VO ) 
X E/ y/ 
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and 


( 12 ) 


• P * 1 /* * \ 


From Eqs. (9) to (12), we obtain, 

E, 


1 


5 -^4 V ^ 




(13) 

(14) 


Now from the given uniaxial loading conditions the strain rate along y is given by 



• _ • e , ♦ p 

y y y 



• e 1 • p 
^x “2 ^x 

(15) 

Substituting 

• p • • e 

e = e - e 
XXX 

(16) 


* * F 

„ a e ^Xo 


and 

• e X _ X 2 

S - E, - E, 

(17) 

into Eq. (15) 

we obtain 



• • r 1 /I N 

^y " ^x 2 ^ 2 ■ 


(18) 


Substituting into Eqs. (7) and (14) and then into Eq. (16), the following thickness 
ratio is obtained, 

E, 


= 5 - 4v 


t^ 4 (1 - v2) ^x 


(^- 1) 


(19) 


From which 


where 


t, + 


, n - 2v)‘ 

5 - 4v 


- Ex2 

^ ^X2 


(20) 


( 21 ) 
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It is noted that this factor 3 is the same for isotropic materials for which 
Eqs. (9) and (10) hold. For one dimensional problem the ratio of the area of 
plastic element to the total area is 

^ “ ^1 ~ ^2 (22) 

A E, 

Thus 3 is the modification factor for plane stress problems. For v = 0.3, of 
Eq. (22) 3 = (l-2v)^/(5 - 4v) = 0.0421. Thus, the modification is very small when 
the uniaxial behavior along the longitudinal direction is used to determine the 
thickness ratio. 

Now if a coupon is cut along the y-axis, then when a = a = o =0, 

X x^ 

a 

^1 

=2a (23) 

^1 

and e, P = - 2a kf (24) 

Yl X-, 


Following the same derivation given above, if the tangent modulus is E^ 
thickness ratio becomes 


2 


the 



^1 = ‘'I ■ ^yz 

‘ - s' 

(25) 

where 

„ ^ _ (1 - 2a v)2 

^ . . z 

(26) 


1 - 4va + 4a 


For V - 0.3 and fora between 0.5 and 2, the values of 3" is in the range of 
0 to 0.35. 


Equating Eqs. (19) and (26) and solving for E-j^ we obtain a relation 


between E and E , 
X 2 y2 



(1 - 3) ^1 ^X2 

(1 -3") E^ - (3 -3") 


(27) 
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Thus for anisotropic material represented by the mechanical sublayer model the 
initial tangent moduli for the two directions are not the same. When the 
stress-strain curves for two perpendicular directions are given, a mechanical 
sublayer can be obtained according to the following procedure. 

(1) For the stress-strain curves, for both x-and y-direction determine the 

initial yield stresses 'i' and 'i' . These are also the yield stresses Y and 

Xi y^ 

Y of the sublayer and are euqal to E, and Et c , respectively. 

y-j I X.J I y-j 

(2) Based on one of the curves say for the x direction try an initial tan- 
gent modulus E . Knowing Y , Y , and E , a and 3 for layer 1 can be 

X2 x^ y-j X2 

determined, hence, t^/t can be calculated from Eq. (2) and the initial tangent 

modulus E for the y-direction, from Eq. (27). 

^2 


(3) The 2-sublayer model is then used to analyze two uniaxial loading 

problems and to obtain the stress-strain curves for both directions. The 

intersections of these lines to the actual stress-strain curves now determines 

the second set of transition points at which the second layer yields. The 

stresses and strains are respectively 'i' and Y"' and e and e as shown 

X 2 Y2 ^2 ^2 

in Figure 3. In general, the yield stresses Y and Y of the new sublayer are 

X 2 y2 

not equal to E, e and E, e , although in the case of subelement model for 
uniaxial problems, the yield stress for the second subelement is equal to 

For 3-D problems with isotropic plastic behavior, Hunsaker has obtained 
closed form solutions for the case of a 2-subelement model. For a typical 
case with E 2 = 0.5 E-j and ~ 3 g-| the difference between Y^ and E-| 
only 5%. 

The choice of the initial tangent modulus E must be such that the second 

X 2 

set of transition points do fall on the actual stress-strain curve. In 
general, an iterative procedure is required. 
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(4) 


The plate is now considered as a new two-layer model with the yield 


stresses for the two directions equal to Y and Y 


^2‘ 


Then by choosing a 


tangent modulus E for the next segment, values of (t-, + tp)/t, E , and Y 

X3 \ c > 

Y" , Y , and Y can be determined following the same steps (2) and (3) 


above using 


^1 h ^1 ~ % 

* 

and E^ = H ~ g) E-, Exg 

^ (1 - e^) E, - (6 - B') E^ 

I X 3 

for the value of B' Eq. (26) is used with a = (Y /Y )^. 

X 2 


(28) 

(29) 


(5) The mechanical sublayer model may be constructed using the stress- 

strain curve for the transverse direction as reference. In that case, Eq. (25) 
as used to determine the thickness ratio's and Eq. (27) is used to solve for 
E in terms of E . 

Xi 

In the solution of a plane stress problem given in this paper a simplified 
procedure was adopted. The stress-strain curves for either the longitudinal or the 
transverse directions is replaced by a curve with linear segments. With one of the 
tangent moduli given, the other tangent modulus can be obtained and is again consid- 


ered constant. Thus, the intercept of this linear segment to the actual stress-strain 
curve can be determined. The yield stresses for Sublayers are obtained by the 
simple formulas, 


and 






^i^y. 

1 


(30) 
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3. 


Fi n i te Element Method Based on Viscoplasticity Theor. 


Time independent elastic-plastic analysis can be made by using viscoplasticity 
models [ref. 7]. In the case of elastic-perfectly-plastic material, it is 
only necessary to consider the rate of viscoplastic strain components as 


•vp ^ A V. 9F 

® = Y < ♦ 


where F (a) = 0 represents the yield surface and 


< c() > 



for a > 
for a = Oy 


(32) 


a = equivalent stress 

0 y = yield stress under uniaxial loading 

Y = a fluidity parameter which is arbitrary for the corresponding 
elastic-plastic analysis. 

A corresponding viscoplastic model for the mechanical sublayer model is an 
arrangement of viscoplastic elements in parallel as shown in Figure 4. 

In the finite element analysis of a plane stress problem, the entire 
panel is discretized into N finite elements, the thickness of each of which is 


divided into M sublayers according to the modelling described in the previous 
section. A finite element method for the creep problem has been formulated by 
the initial strain approach using the assumed stress hybrid model [8,9]. 

The method is extended to the present multilayer model. For the incremental 
solution of the elastic-plastic problem, the procedure is as follows, 

(1) An elastic solution of nodal displacements due to a given load 
increment is made using the assumed stress hybrid elements. 
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(2) The stress increments at selected Gaussian stations in each 
element are evaluated, 

(3) The stresses in each sublayer at each Gaussian station are updated 
and the corresponding equivalent stress o’ evaluated, 

(4) A time increment A t is chosen and the increments of viscoplastic 
strains are Ae'^^ = A t for all sublayers at all Gaussian stations. Here the 
viscoplastic strain rates are determined by Eqs. (31) and (32). 

(5) The equivalent nodal forces due to the viscoplastic strain increments 
can be evaluated. They are used for the determination of nodal displacements for 

the time increment. 

(6) Steps (2) and (5) are repeated until the changes in stresses within the 
time increment become less than a small prescribed limit. 

At this stage, the equivalent stress o in each sublayer is either equal or smaller 
than the yield stress of that sublayer and stabilized stress state for this 
loading increment is obtained. The use of successive time increments is 
equivalent to an iterative procedure for this elastic-plastic problem. 

A guideline for the choice of the time increment At to assure numerical 
stability is, according to Cormeau [ref. 10], 

At £ 4 (1 + v) 

3 E Y (33) 

where E is the Young's modulus 

4. Anisotropic Analysis of Shear Lag Structure 

To illustrate the finite element solution using the present anisotropic 
model, a shear lag structure which was tested in 1963 at Massachusetts Institute 
of Technology is used [ref. 11]. The structure shown in Figure 5 was integrally 
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machined from thick 2024-T4 aluminum alloy plate. It is a rectangular panel 27.94 cm 
X 25.4 cm X .203 cm, stiffened by tapered stiffners along the loading (y) 
axis. The stiffness became flush with the panel at the center, where the 
stresses and strains are the highest under the applied tension load through the 
stiffners . 

Fitted stress-strain curves for both the longitudinal or x-direction and 

the transverse or y-direction were obtained from tension tests as shown in 

Figure 6. An average of these two curves is also shown. 

A mechanical sublayer model for this material was constructed by the 

simplified procedure indicated in section 2. The corresponding curves with 

piecewise linear senments are shown in Figure 7. The yield stress under 

shear was not obtained in the original experiment. For the mechanical sublayer 

model, the values of Y for each layer is set equal to Y/ /3 where Y is the 

xy 

larger of Y and Y . Of the five sublayers used in the models the last one is 
X y 

considered elastic. The strain distribution of the panel, thus, is determined 
directly by the elastic strains of this layer. Table 1 lists the thickness ratios 
and yield stresses of the various sublayers. 

For the finite element analysis one quarter of the panel is used. It is 
subdivided into 7x7 rectangular plane stress elements as shown in Figure 8. 

The tapered stiffness is also modelled as plane stress elements with constant 
thickness in each element. For the finite element solution all numerical 
integrations were obtained by using 2x2 Gaussian quadrature. 

The resulting strains and at the center of the panel were determined 
by extrapolating from the two Gaussian stations along the diagonal line of 
the element at closest to the center. Their variations with respect to the 
applied load are shown in Figure 9. Also plotted for comparison are: 
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(1) A finite element solution obtained by the present method but by mod- 
elling as an isotropic material based on the average stress-strain curve shown 
in Figure 6. 

(2) A finite element solution obtained by Jensen et al . [ref- 12], using 

144 constant strains triangular elements and by modelling as an isotropic material. 

(3) Experimental results obtained in Ref. 11. 

It is seen that the solutions by the three finite analyses are comparable 
although they do not agree with the experimental result, especially, for 
the component. In this particular case the approximate solution obtained 
by considering the material property to isotropic appears to be very close to 
that by the present modelling of anisotropic material. 

5. Conclusions 

A method has been developed for the modelling of anisotropic plastic 
behaviors for metals by the mechanical subelement model. It has been incor- 
porated in a finite element analysis program based on the assumed stress hybrid 
model and, on the viscoplasticity theory. 
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Table 1. 


Constants for Hechanical Sublayer Model (Anisotropic Plasticity) 













































TENSILE STRAIN 10' 


Figure 6. Fitted Stress-Strain Curves of 2024-T4 Aluminum Alloy Plate 



Figure 7. Approximate Stress-Strain Curves by Mechanical Sublayer Model (2024-T4) 



-H-h- 


i-+-+ — 

I I I 

't“t i* “H ^ 

^ 


T 

2a 

i 

4 

2a 

4 - 




/// 7 y ' //// / o.i 

* • a 4 a* *' a 2a 2a 3a » ■»■ ^~ 


I 


6a 


a = 1.27 cn 

Figure 8. Finite Element Mesh for Shear Lag Structure 
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Figure 9. Strain Components at Center of Shear Lag Structure - Comparison of 
Different Finite Element Solutions 
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